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THEORY OF LOW TRANSITIONS IN CO DISCHARGE LASERS 


Barry D. Sidney, Roy M. Mclnville, * Neill S. Smith, * 
and H. A. Hassan* 

Langley Research Center 

SUMMARY 

A self-consistent theoretical model which couples the electron and heavy particle 
finite -rate kinetics with the optical and fluid dynamic processes has been employed to 
identify the various parameters and explain the mechanism responsible for producing 
low-lying transitions in slow-flowing CO lasers. It is found that lasing on low-lying 
transitions can be achieved at low temperatures for low pressures (or low -flow rates) 
together with high partial pressures of He and N 2 . The role of N 2 has been identified as 
an additive responsible for reducing the electron temperature to a range where the trans- 
fer of electrical power to the lower vibrational modes of CO is optimum. 

INTRODUCTION 

Molecular gas lasers have received considerable attention in recent years because 
of their high efficiencies and high power capabilities in the infrared region. The CO laser 
is of particular interest because of its high quantum efficiency and its ability to oscillate 
on a number of vibrational -rotational bands. Transitions involving higher vibrational 
bands are strongly attenuated by atmospheric water vapor; however, low-lying transitions 
with wavelengths below 5 jUm propagate in the atmosphere with relatively high transmit - 
tivity. (See ref. 1.) Moreover, the 1 — 0 transition can be used for detecting CO in 
long-path absorption measurements. Thus, the ability of the CO laser to operate on low- 
lying vibrational bands in the 5 -pm region makes it an indispensable tool in atmospheric 
transmission and probing experiments. 

The objective of this study is to identify the various parameters and explain the 
mechanism responsible for producing low-lying vibrational transitions in CO laser 
systems. Experiments (refs. 2 to 4) have shown that low temperature and low CO con- 
centration in a CO-N 2 -He electric discharge system are two factors influencing the 
appearance of low-lying vibrational transitions. Both of these effects were evident in 
the recent analysis of the high -flow laser systems by Smith and Hassan. (See ref. 5.) 


* 


North Carolina State University at Raleigh. 


For the study reported in this paper, a self-consistent model was formulated and was 
used to evaluate the manner in which the operating conditions, that is, pressure, flow 
rate, mixture ratios and power, affect the low-lying transitions. The need for a self- 
consistent model became evident because earlier analyses of slow-flow CO lasers 
(refs. 2 and 6 to 8) assumed certain flow properties which could be determined as part of 
the solution. Thus, the theories of references 2, 6, and 7 assume the electron distribu- 
tion function to be Maxwellian at some given electron temperature; moreover, these 
theories assume the electron nvimber density and the gas temperature. The effect of 
such assumptions is that important phenomena resulting from the coupling of the kinetics 
and fluid dynamic processes of the electrons and heavy particles are ignored. Center 
and Caledonia (ref. 8) have also presented calculations for slow-flow CO lasers in which 
the distribution function was obtained from a solution of an appropriate Boltzmann equa- 
tion; however, they had to assume a value for the gas temperature. In contrast to these 
theories, the present model couples the kinetics of the electrons and heavy particles with 
the optical and fluid dynamic processes in the laser system. It treats molecules in dif- 
ferent quantum states as different species and employs the electron Boltzmann equation 
and the conservation of species, mass, momentum, and energy equations to describe the 
system. For a given geometry, voltage, current, gas composition, pressure, mass flow 
rate, and reflectivities, the theory can predict the number densities of the electrons and 
the excited states of the heavy particles, the electron and gas temperatures, the intensi- 
ties of the lasing transitions, and efficiency. 

These theoretical results are compared with gain and intensity measurements 
obtained with a slow-flow CO electric discharge laser of conventional design. (See ref. 9 
and fig. 1.) In general, the comparisons show good agreement and indicate the validity of 
the model developed here. 

The work carried out at North Carolina State University at Raleigh was supported in 
part by NASA Grant NSG-1013. 


SYMBOLS 

A cross-sectional area of laser tube 

Agt empirical constant in equation (A20) 

A(u) quantity defined by equations (B7) 

Ay^^A^^^ spontaneous emission coefficients of vibration level v of CO 
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a radius of laser tube 

aQ Bohr radius 

BgQ rotational constant for ground state of species s 

B{u) quantity defined by equations (B7) 

BVTy backward vibration-translation rate coefficient defined by equation (A17) 
By rotational constant of vibrational level v of CO defined by equation (C8) 

bgj. quantity defined in equations (A3) 

C(u) quantity defined in equations (B7) 

c velocity of light 

D(u) quantity defined in equations (B7) 

Dy quantity defined in equation (C9) 

dg,d^ collision diameter of Lennard-Jones potential for species s and t 

dg|. collision diameter defined in equation (A8) 

E electric field strength 

E electric field strength vector 

s t 

E^jEy energy of vibrational level v of species s and t 
e electronic charge of electron 

F adiabaticity factor defined by equations (A4) and (A5) 

Fy quantity defined by equation (C7) 

fg collision frequency defined in equations (C4) 



electron energy distribution fxmction 


f(u) 

Gs 

g(u) 

H 

Ho 

h 

*^s,v 

I 


Is.It 

J 


j 

I 


Kd 


Kq 

KsD 


K 


K 


K 


sQ 

st 

VT 

st 

VV 


anharmonic vibrational energy function defined by equations (16) 
function defined by equation (Bll) 

specific enthalpy of gas mixture defined by equation (15) 
specific enthalpy of gas mixture at tube entrance 
Planck's constant 
specific enthalpy of species s 

specific enthalpy of vibrational level v of species s 
arc current 

ionizational potential of species s and t 
intensity 

rotation quantum number 
current density 
current density vector 
quantity defined in equations (B8) 
quantity defined in equations (B8) 
quantity defined in equations (B3) 
quantity defined in equations (B5) 
empirical constant in equation (A2) 
empirical constant in equation (A21) 
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k 


Boltzmann's constant 


L length of laser tube 

Lgpt length of optical axis 

, long-range vibration -vibration rate coefficient defined by equation (A20) 

v,v 

empirical constant in equations (A3) 

M average mass defined in equations (B8) 

mg electron mass 

mg mass of species s 

m^. quantity defined by equations (A3) 

rh total mass flow rate 

riig mass flow rate of species s 

N total number density 

Ng density of electron 

Ng number density of species s 

Ng j number density of species s with rotational quantum number J 

Ng y number density of species s in vibrational level v 

P(J) rotational branch of a given vibrational transition 

Py power from vibrational level v 

P® T rate for electron impact excitation from vibrational level J 

V.d 
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p 

P^st 

Qs 

Qsj 

^sm 

Qsv 

Q»P 

Qw 

QpQ2 

q 

RgjRt 

^s,v 


Sv 


electron inapact vibrational excitation rate for species s defined by 
equations (B20) and (B21) 

static pressure 

relaxation time of species s by species t 

energy loss due to lasing 

energy loss by spontaneous emission 

electron inelastic cross section with energy loss Ugj for species s 

electron inelastic cross section with energy gain Ugj for species s 

momentum transfer cross section for electrons in species s 

electron inelastic cross section with energy loss Ug^ for species s 

total heat loss 

heat loss by wall transfer 

cross sections defined in equations (B8) 

heat flux vector 

particle rate of production of species s and t 

particle rate of production of vibrational level v of species s 

vibrational matrix element of species s 

ionizational energy of ground state of Bohr atom 

distance in radial direction 

mirror reflectivities 

line shape factor defined by equation (C3) 
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vv; 


st 
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st 


v,v' 


V>J1 


short-range vibration -vibration rate coefficient defined by equation (A21) 

species 

temperature 

critical temperature used in equation (AlO) 
electron temperature 

characteristic temperature defined in equations (A3) 

vibrational temperature of species s defined by equation (8) 

wall temperature of laser tube 

flow velocity 

vector flow velocity 

mean flow velocity 

electron energy in electron volts 

electron energy loss or gain in Jth inelastic collision process with species 
gas temperature in electron volts defined in equations (B8) 
arc voltage 

vibration-translation rate coefficient defined by equations (A2) and (AlO) 
vibration-vibration rate coefficient defined by equation (A19) 
vibration quantum number 

vibration distribution truncation level of species s 
number of lasing transitions 


Illlllllllllllll 


Wg^ empirical constant used in equation (A20) 

w electron drift velocity defined by equation (B15) 

Xg mole fraction of species s 

X quantity defined in equations (C4) 

pNg 

Yg mass fraction of species s, 

Yg y mass fraction of species s in vibration level v 

y variable used in equations (A4) and (A5) 

y^ quantity defined in equations (A3) 

st 

y^ quantity defined in equation (A22) 

Zg^ collision frequency defined by equation (A7) 

0 !^ threshold gain for transition v 

/3g mole fraction defined in equations (B8) 

r gamma function 

Yy j gain coefficient defined by equation (Cl) 

AE energy loss found in equations (Al) and (A18) 

At'g Lorentz half -width defined in equations (C4) 

Aiq^ Doppler half -width defined in equations (C4) 

5g quantity defined in equations (A3) 

eg,6j. potential well depth of Lennard- Jones potential of species s or t 

€g^ potential well depth defined in equation (A6) 
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.St 

®v,v* 


variable used in equation (A23) 
electric dipole moment of species s 
electric quadrupole moment of species s 
characteristic temperature defined in equations (A3) 
temperature defined in equation (A22) 
quantity defined in equations (B8) 

mean thermal conductivity of mixture defined by equation (18) 

thermal conductivity of species s 

electron mobility defined by equation (B16) 

reduced mass defined in equations (A3) 

frequency defined by equation (C2) 

frequency defined by equation (C6) 

density 

optical cross section 

quantity defined in equations (B3) 


CTgQ quantity defined in equations (B5) 

(2 2 ) 

Or ’ ’ quantity defined by equation (A9) 

) vibrational constants in equations (16) 

An arrow over a symbol denotes a vector quantity. An asterisk denotes the largest 
possible value. A prime denotes a variation of any term. The subscript max denotes 
a maximum value. 
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ANALYTICAL FORMULATION 


The laser system modeled in this analysis is based on the experimental work 
of Sidney et al. (ref. 9) and is shown schematically in figvire 1. The discharge tube is 
made of pyrex and is 2.54 cm in diameter and' 92 cm long. Annular stainless-steel 
electrodes, 2.54 cm in diameter and 8 cm long, are fitted to the tube ends and calcium 
fluoride Brewster angle windows are placed on the electrode ends. The tube, exclu- 
sive of the electrodes, is immersed in a liquid nitrogen bath in a bakelite container. 
Compartments at either end of the container were designed with adjustable liquid nitro- 
gen leaks to control the temperature of the electrodes. The constituent gases are pre- 
mixed in an upstream reservoir and are introduced at the cathode and exhausted near 
the anode. The pressure drop in the tube was negligibly small for the operating con- 
ditions of the experiment. 

In analyzing the experiment, molecules in different vibrational quantum states are 
treated as separate species. For such a multicomponent mixture, the governing equa- 
tions are the electron Boltzmann equation, the conservation of species equations, the 
overall continuity, momentum, and energy equations, and the equation of state. 

The electron Boltzmann equation is 


e 2 d 

u df 

A 

■* /u2n Q a 

3 du 

^(NsQsm) 

_s 

^ L 

s 

mg dul ^s^sm^j 



!!!lekT^/ 2n o 

me e dur ^s^sm 


s 


Ng/^ [(u + Ugj) f(u + Ugj) Qgj(u + Ugj) 


- U f(u) ^ [Qsj(u)] + J [(u - Ugj) f(u - Ugj) Q*j(u - Ugj) 
J J 



This equation is considered in detail in appendix B. The electron temperature and the 
electron number density are determined from the current density and the electron dis- 
tribution function as 
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( 2 ) 


T 


e 



u 


3/2 


f(u) du 


and 


N 


e 


_L 

ew 


w= mE = mn| 



where 


M = 




(3) 


By use of the steady -state approximation, the conservation of species equation for 
a given species s reduces to 


Rs = 0 


(4) 


The species considered are the vibrational levels of CO from v = 0 to 
the vibrational levels of N 2 from v = 0 to v = where v^q and 
chosen to be 50. 


V = ^CO 
Vn 2 were 


The kinetic processes assumed to be of importance in developing expressions for 
the production rates Rg are 

(1) Electron impact excitation of vibrational levels 1 to 8 of both CO and N 2 

(2) Single quantum vibration -vibration (V-V) exchange collisions 

(CO-CO, CO-N 2 , N 2 -N 2 ) 

(3) Single quantum vibration -translation (V-T) exchange collisions 

(CO-CO, CO-He, CO-N 2 , N 2 -CO, N 2 -N 2 , N 2 -He) 

(4) Spontaneous emission of CO (Av = 1 and Av = 2) 

(5) Stimulated emission and absorption of CO. 
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The cross sections for electron excitation of vibrational levels 9 and 10 for both CO 
and N 2 were measured by Boness and Schulz. (See ref. 10.) However, because of their 
small magnitudes, they were not included in this analysis. 


For the vth vibrational level of CO, 


^CO,v " ^e 


8 

8 


1 (nco,j4°) 

1 

0 

p 
< 0 

c-iO 

^CO,v+l 

J=0 

J=0 

-1 


CO 
j=o 




J=0 




- N, 


CO,v 


^CO 


Z ('•co. 


J=0 


* 

VN. 


2 

/ No-CO UU-JNo\ 

^ Z N j-^Z'^spVT— + VT, 

J=0 


,„,.CO-CO , T,T iM7^CO-CO\ 


CO-N. 


CO-s . ,,^rr>CO-sj 




* 

VNc 


+ N, 


CO,v-l 


CO 

Z rCO,J+l''''j+l,v-V 

J=0 


J=0 S 


+ o A^2) I yy+l,J^+l,J _ 

CO,v+2 v+2 hi^v+l,J j 


(5) 
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and for the vth vibrational level of N£ 

8 , „ , 8 




t (%2. a') - %.v 1 ('’vj) 

J=0 J=0 


+ N. 


N 2 ,v +1 


* 


Z V^N 2 ,J v+l,j/ 


★ 

vco 


J=0 S 


N, 


N 2 ,v 


J=0 




*2 


I K2, 


N 9 -N 9 


♦ 

^co 


N 9 -N 9 \ V ( CO-No No-CO\ 


J=0 


I N 9 -S N 9 -S 

^Z'^sl^VT^+y ■^VT/ 

s 


+ N. 


N 2 ,v -1 


N 2 

sK 

J=0 ' 


N 9 -N 9 

j+ivvj^y,v!i 


* 

^CO / 


t (nco.j^iVv“i;v-ZE(’‘3®''t”= 

J=0 s 


-s 


(6) 


These expressions contain terms which involve N^q 1 ^ 

terms are simply set to zero, then false inversions can occur on the upper levels. To 
eliminate this problem, the number densities for the v^q + 1 and the + 1 levels 

are estimated by assuming that the population distribution in the last few levels can be 
approximated by a Boltzmann distribution; that is. 


Ns,v*+1 = Ng^v* 


E®+ . - E®*) 

Vg+1 Vg 


kT 


(7) 


s,v 


where 


, -(<: - 

■s,v 

" '\<-y 


( 8 ) 


k log, 


This method was first used by Lordi et al. (See ref. 11.) 
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The various V-T and V-V rate coefficients used in equations (5) and (6) are dis- 

g 

cussed in appendix A. The electron vibrational excitation rates Pj^j are defined in 
appendix B. The only excitation cross sections which have been measured experimen- 
tally for CO and N 2 are those for excitation. from the ground state (ref. 10). Chen 
(ref. 12) has calcvilated theoretical cross sections for vibrational excitation from levels 
V = 1, 2, 3, and 4 for N£. By using his results, Abraham and Fisher (ref. 13) have dem- 
onstrated that if electron vibrational excitation processes have the same Av, then they 
have approximately the same rate as the excitation process from v = 0. Thus, it is 
assumed that 


V, v+1 



(v = 1 to 7) (9) 


The spontaneous emission coefficients for CO were generated from the matrix ele- 
ments of Young and Eachus (ref. 14). However, they calculated the matrix elements only 
up to V = 24. The elements for the higher levels were generated by extending the calcu- 
lations of reference 14 to v = 50. The dipole moment vibrational elements and the 
spontaneous emission coefficients are given in table I. 

The overall continuity equation is written as 


pUA = m 


( 10 ) 


whereas the overall momentum equation reduces to 
VP = 0 


( 11 ) 


The overall energy equation has the form 


V 


pUH 




E 


( 12 ) 


For the steady -state approximation employed here, this equation can be integrated to 
give 

m(H - Hq) = IV - Qt (13) 


The total heat loss Qj- is expressed as the sum of losses due to conduction to the wall, 
spontaneous emission, and lasing. Thus, equation (13) becomes 

m(H - Ho) = IV - - Qs - Qf (14) 
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TABLE I, - DIPOLE MATRIX ELEMENTS AND SPONTANEOUS 
EMISSION COEFFICIENTS OF CO 


V 

rCO /„C0 2 
^v,v-l/^l,0 

'^-1 

sec ^ 

a(2) 
sec ^ 

V 

rCO Leo 2 
^v,v-l/^l,0 

^-1 

sec 

A2) 

^ -i 

sec ^ 

1 

1.000 

33.4 


26 

19.698 

227.6 

201.8 

2 

1.995 

64.2 

0.9 

27 

19.987 

220.0 

212.2 

3 

2.985 

92.5 

2.7 

28 

20.241 

211.9 

222.2 

4 

3.968 

118.4 

5.3 

29 

20.460 

203.7 

231.9 

5 

4.945 

141.9 

8.8 

30 

20.644 

195.4 

241.1 

6 

5.913 

163.3 

13.0 

31 

20.793 

186.9 

249.9 

7 

6.867 

182.3 

18.0 

32 

20.908 

178.4 

258.2 

8 

7.816 

199.4 

23.7 

33 

20.988 

169.8 

265.9 

9 

8.739 

214.2 

30.2 

34 

21.033 

161.3 

273.0 

10 

9.660 

227.4 

37.2 

35 

21.044 

152.9 

279.6 

11 

10.560 

238.7 

45.0 

36 

21.020 

144.5 

285.5 

12 

11.436 

248.0 

53.4 

37 

20.961 

136.3 

290.7 

13 

12.285 

255.5 

62.3 

38 

20.867 

128.2 

295.3 

14 

13.104 

261.3 

71.6 

39 

20.738 

120.4 

299.2 

15 

13.898 

265.5 

81.4 

40 

20.575 

112.7 

302.4 

16 

14.656 

268.2 

91.6 

41 

20.377 

105.2 

304.9 

17 

15.368 

269.2 

102.3 

42 

20.144 

98.0 

306.8 

18 

16.038 

268.8 

113.0 

43 

19.877 

91.0 

307.9 

19 

16,644 

266.8 

124.1 

44 

19.575 

84.2 

308.3 

20 

17.220 

263.8 

135.3 

45 

19.238 

77.8 

308.0 

21 

17.735 

259.6 

146.6 

46 

18.866 

71.6 

307.1 

22 

18.227 

254.8 

157.8 

47 

18.459 

65.6 

305.5 

23 

18.676 

249.1 

169.1 

48 

18.017 

60.0 

303.2 

24 

19.056 

242.4 

180.2 

49 

17.540 

54.7 

300.4 

25 

19.375 

235.0 

i 

190.8 

50 

17.028 

49.6 

296.9 



The specific enthalpy of the mixture is given by 


XO 


* 


H * Wne + 


v=0 


Yxt _.h- 


CO,v“CO, 


v=0 


N2,v“N2,v 


(15) 


where 


h - 5 ^'^e , ^CO 
® 2 me nie 

5 kT 


‘He “2 


«iHe 


7 kT ^ 


CO 


CO,v 2 mcQ m^o 

E^2 

, _ 7 kT ■*^v 

%,v- 2 


(16) 


mN2 niNg 

K = - Gs(0)] 

Gs(v) = 0Je(v + I) - WeXe(v + + «eye(v + + ‘^eZe(v + 1)^ 


The units of Gg(v) are cm“^ and the vibrational constants for CO and N 2 foimd in equa- 
tions (16) are given in table II. 

TABLE II. - VIBRATIONAL CONSTANTS OF CO AND N 2 
[From refs. 15 and 16] 



CO 

N2 

"e 

2169.82 

2359.61 


13.292 

14.456 

"eYe 

.01082 

.00751 


.000052 

-.000509 
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Heat loss to the wall can be expressed as (ref. 17) 

Qw “ 67 t(T - T^)XjnL 


(17) 



where 


X 


m 






(18) 


Equation (17) is developed in reference 17 on the assumption that the average or bulk 
temperature of a fully developed flow in a tube of radius a is defined as 


T = I- \ T(r) dr 
a Jq 

But the usual definition of the average gas temperature (ref. 18) is 


(19) 


J T(r) U(r) dA 
J U(r) dA 


( 20 ) 


Thus, the assumption of a constant U(r), which is implied in reference 17, would seem to 
indicate that this expression for is incorrect. However, for fully developed pipe 

flow (ref. 18), 


U(r) = Uni 




( 21 ) 


Using equation (21) in equation (20) leads to the expression for in equation (17). 

The thermal conductivities used in equation (17) were obtained from reference 19. 


Spontaneous emission losses are given by 


Qs = AL 


* 

^CO 


* 

VCO 


I NcO,v 4*’K° - <-°l) * I - EvS) 


V=1 


v=2 


( 22 ) 


and the laser power output is 

«« = AL J J I, J (23) 

Finally, the equation of state for the system can be written as 

p = kNgTg + kNT (24) 


17 



When operating in the lasing mode, the intensities of each transition 1^^ j assume 
values such that the gain coefficient j for each transition satisfies the threshold 
condition (ref. 20) 


-logg (rjrg) 

~ 2L„pt 


(25) 


This condition is a statement of the fact that the intensity gain on a line just equals the 
intensity loss on that line for a round trip through the discharge region. If the losses 
are of a magnitude such that the gain coefficient for a given transition is less than the 
threshold value, no lasing can occur on that line. 

An expression for the gain coefficient including the effects of both Doppler and 
pressure broadening is given in appendix C. The development of this expression includes 
the assumption of the equilibrium of the rotational states at the gas temperature. Thus, 
it is not possible to predict transitions on more than one rotational line in a given vibra- 
tional band. Therefore, it is assumed that for a given vibrational band, lasing occurs 
only on the P(J) branch corresponding to the maximum small signal gain coefficient. 

Equations (1), (4), (10), (11), (14), and (24) together with the expressions for the 
electron temperattire and the electron number density given by equations (2) and (3) form 
the basic system to be solved. When the system operates in an oscillator mode, this 
system of equations is supplemented by threshold conditions given by equation (25) equal 
in number to the lasing transitions. 


METHOD OF SOLUTION 


The governing system of equations consists of the electron Boltzmann equation, 
which is a difference differential equation, and 3 + Vqq + Vj^^ + algebraic 

equations corresponding to the number densities of the Vqq vibrational levels of CO 
plus the ground state, the number densities of the vj^g vibrational levels of N£ plus 
the ground state, the gas temperature, and the v-p lasing intensities. For small signal 
gain calculations, the intensities are set equal to zero and thus the last vp equations 
are not needed. The electron Boltzmann equation is integrated numerically by using a 
constant -step -size, fourth-order, Adams-Moulton predictor -corrector method. The sys- 
tem of nonlinear algebraic equations is solved numerically by using a Newton -Raphson 
technique. 

The procedure used in obtaining a complete solution for a given set of operating 
conditions is summarized as follows: 

(1) For the given composition, the electron Boltzmann equation is solved for a 
range of E/N values to obtain the electron distribution function. This solution is then 
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used to calculate the electron vibrational excitation rates which are needed in the rate 
equations. The electron temperature and mobility are also calculated. This information 
is stored for later use. 

(2) Next, a calculation of the small signal gain coefficients is carried out. With the 
given operating conditions, that is, composition, pressure, flow rates, arc voltage and 
current, and tube dimensions, an initial value for the temperature and an initial vibra- 
tional distribution for the CO and N 2 molecules are assumed. The corresponding electron 
temperature is determined by interpolation from the table generated in step (1). A simi- 
lar procedure is employed to determine the vibrational excitation rates and mobility. 

The electron ntimber density can then be found from equation (3). The V-V and V-T 
rates for the various collisions are computed at the assumed temperature. By using this 
iitformation, the initial values of the rate equations and the energy equation together with 
all their partial derivatives which are needed for the Newton -Raphson iteration scheme 
are calculated. After each iteration, the procedure of determining an electron tempera- 
ture and all parameters associated with it and a set of V-V and V-T rates for the new 
temperature is repeated. The iteration is continued until convergence on all variables 
to a relative error of 10"^ is achieved. 

For the small signal gain calculations, all intensities are set to zero and thus 
the equations for the threshold condition are omitted. The converged values of number 
densities and temperature are used to calculate small signal gain coefficients for the 
vibrational -rotational transitions of interest. For each vibrational band, the gains are 
calculated for a number of P(J) transitions and the value of the maximum gain on each 
band together with the J number of the rotational transition on which it occurs are 
retained for use in step (3). 

(3) Finally, an intensity calculation is carried out. In this case, it is necessary to 
introduce the threshold condition given by equation (25) for each intensity to be calculated. 
Thus, v-j. additional unknowns and v-j. additional equations are added to the system. 

The lines for which small signal gain coefficients are maximum for a given band 
may oscillate if their gain coefficients are greater than or equal to the threshold value. 
Such lines are determined by comparing the maximum gain values for each band calcu- 
lated in step (2) to the threshold gain value. 

The resulting system of equations is solved again by a Newton -Raphson technique. 

If a negative intensity appears in the converged solution, it is dropped, together with its 
associated threshold condition, and the solution procedure is repeated. When a con- 
verged set of all positive intensities is obtained, the power from each line is calculated 
from the expression 
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The power from all lines is summed to obtain the total laser power. 

An important point to remember when calculating intensities is that the J num- 
ber which specifies the P(J) transition with the maximum small signal gain coefficient 
on a given vibrational band is not necessarily the same as for the case when lasing is not 
present. Thus, after each iteration, a small signal gain calculation is performed and 
the J numbers associated with each lasing transition are adjusted accordingly. A typi- 
cal example of this shift can be seen in table HI corresponding to the conditions for 
case 1. (See table IV in "Results and Discussion.") 

TABLE m.- COMPARISON OF J VALUES FOR SMALL 
SIGNAL GAIN CALCULATIONS AND INTENSITY 
CALCULATIONS FOR CASE 1 


Vibrational 

band 

Gain 

calculations 

Intensity 

calculations 

4-3 

P(10) 

P(10) 

5-4 

P(8) 

P(10) 

6-5 

P(7) 

P(10) 

7-6 

P(7) 

P(9) 

8-7 

P(7) 

P(9) 

9-8 

P(7) 

P(8) 

10-9 

P(7) 

P(8) 

11 - 10 

P(7) 

P(8) 

12 - 11 

P(7) 

P(8) 

13 - 12 

P(7) 

P(8) 

14 - 13 

P(7) 

P(7) 


As was indicated earlier, 50 vibrational levels were carried for both CO and N 2 in 
all calculations. Provision was made for calculating a maximum of 25 lasing intensities. 
However, typically, it was not necessary to calculate more than 14 corresponding to 
vibrational transitions from the 3 — 2 up to the 16 — 15. 

RESULTS AND DISCUSSION 

The primary processes of pumping energy into the CO vibrational levels are through 
electron excitation of the lower levels and V-V pvunping of the higher levels. Radiative 
decay and V-T collisions are the primary processes for deactivating the vibrational 



states, and these processes are dominant for the higher levels. Thus, to promote lower 
level transitions, one needs to increase the efficiency of the lower level excitations and 
V-T collisions and reduce the efficiency of V-V pumping. To reduce the efficiency of 
V-V pumping, one needs to reduce the pressure (or flow rate) and partial presstire of CO 
in the system. The presence of a large fraction of He increases the efficiency of the 
V-T collisions and helps maintain a low gas temperature. Nitrogen was found to be bene- 
ficial in promoting low level transitions in the experiments of Djeu (refs. 3 and 4); he 
attributed its role to the fact that v — v - 1 transitions in CO are in resonance with 
the V + 6 — V + 5 transitions in N 2 . 

All calculations performed in this analysis are based on a wall temperatvure of 77 K 
and mirror reflectivities of 1.0 and 0.84. The other parameters employed in the calcu- 
lations, which correspond to the range of operating conditions employed in the experi- 
ment (ref. 9), are listed in table IV. 

TABLE IV. - OPERATING CONDITIONS FOR THE CASES CONSIDERED 


Case 

number 

10"^ g/sec 

10"^ g/sec 

mHe> 

10"'^ g/sec 

Arc 

current, 

mA 

Arc 

voltage, 

kV 

Tube 

pressure, 

torr 

1 

1.96 

0.00 

3.91 

7.3 

3.4 

.3.5 

2 

5.85 

.00 

11.73 

6.6 

4.8 

6.3 

3 

9.77 

.00 

19.55 

6.0 

6.0 

8.2 

4 

1.92 

3.81 

3.83 

2.0 

4.0 

3.5 

5 

5.76 

11.42 

11.50 

6.4 

5.2 

6.3 

6 

9.60 

19.04 

19.16 

5.7 

6.6 

8.2 

7 

1.68 

26.78 

3.38 

6.3 

5.4 

3.5 

8 

5.04 

80.36 

10.15 

5.3 

7.4 

6.3 

9 

8.39 

133.90 

16.92 

4.4 

9.2 

8.2 

10 

1.68 

26.80 

3.88 

6.6 

4.8 

3.8 

11 

5.04 

80.40 

10.20 

5.4 

7.2 

6.8 

12 

8.39 

134.00 

16.90 

4.2 

9.6 

9,0 

13 

3.32 

26.60 

3.35 

6.5 

5.0 

3.8 

14 

9.96 

79.60 

10.00 

5.2 

7.6 

6.8 

15 

16.60 

133.00 

, 16.80 

3.9 

10.2 

9.0 

16 

6.54 

26.20 

3.31 

6.4 

5.2 

3.8 

17 

19.60 

78.80 

9.93 

5.1 

7.8 

6.8 

18 

32.70 

131.00 

16.50 

3.8 

10.4 

9.0 

19 

12.60 

25.30 

1 3.19 

6.4 

5.2 

3.8 

20 

37.90 

75.80 

1 9.57 

5.1 

7,8 

6.8 

21 

63.20 

126.00 

16.00 

3.7 

10.6 

9.0 

22 

23.80 

23.80 

3.00 

6.1 

5.8 

3.8 

23 

71.20 

! 71.20 

8.98 

4.9 

8.2 

6.8 

24 

119.00 

119.00 

15.00 

3.2 

11.6 

9.0 


21 
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The effect of N2 fraction on the maximum small signal gain coefficient is seen in 
figure 2. The numbers associated with the various points indicate the P-branch transi- 
tion at which the gain coefficient is maximum. As is seen from figure 2, N2 has the over- 
all effect of increasing the efficiency of all vibrational excitations. The primary reason 
for this behavior is attributed to the decrease in the electron temperature from 28 000 K 
in the absence of N2 to 13 600 K in the presence of N2 as is seen in figure 21. This 
decrease in the electron temperature results in more electrical power being transferred 
to the lower vibrational levels of both CO and N2. (See ref. 21.) In addition, as the num- 
ber of N2 molecules is increased, the calcvilations indicate that the percentage of excited 
N2 molecules for v > 3 increases; this increase helps pump the lower CO levels because 
V + 6 — V + 5 transitions in N2 are in resonance with v — v - 1 transitions in CO. This 
mechanism, however, is secondary to the effect of the electron temperature. 

As is seen in figure 2, increasing the N2 fraction favors the appearance of transi- 
tions at low-lying bands and low rotational numbers. This condition explains the behavior 
of the relative intensity indicated in figure 3 which shows low-lying transitions are favored 
with increasing N2 concentration. 

The effect of CO on low-lying transitions is discussed next. Increased CO fraction 
is accompanied by an increase in the V-V pumping or energy transfer from lower to 
upper levels. This increase dominates all other mechanisms responsible for pumping 
the lower levels and is the reason for the decrease in the maximum gain coefficients of 
the low-lying bands with increase in CO concentration shown in figure 4 and the disap- 
pearance of the low-lying transitions in figure 5. 

Figures 6 and 7 show the effect of pressure (or, equivalently, the flow rate) on the 
maximum small signal gain coefficients and the relative intensities. An increase in the 
pressure for a given mixture results in an increase in the number of CO molecules in the 
system and thus improves the efficiency of the CO V-V pumping mechanism. The calcu- 
lations also indicate that at higher pressures, more electrical power is being transferred 
to the lower levels of CO while at the same time the percentage of excited N2 molecules 
decreases. Thus, N2-CO V-V pumping is reduced. These considerations are responsible 
for the decrease in the maximum small signal coefficient and the relative intensities of 
the low-lying levels as the pressure increases. 

In figures 8 and 9 a comparison of theory with experiment is presented. Figure 8 
shows a typical plot of the relative intensity as a fimction of the upper vibrational levels; 
the numbers associated with the various points indicate the P-branch at which the inten- 
sity is maximum. In general, the predictions of the levels at which the peak intensity 
occurs, the relative magnitudes of the intensities, and the rotational numbers at which 
the intensity is maximum in a given band are in good agreement with experiment. 
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However, comparison of measured and calculated gain coefficients for the low-lying 
line 5—4(9) indicated in figure 9 shows that the theory overpredicts the measurements by 
a factor of two. The discrepancy between calculated and measured gain coefficients in 
CO lasers was discussed by Greene and Harris. (See ref. 22.) They showed that their 
calculations of the small signal gain coefficient for low temperature e-beam sustained 
continuous wave CO-Ar lasers, using the expression given in reference 7, agree with 
experiment if they discard hard sphere models for the optical broadening cross sections 
and employ instead cross sections which have the dependence on the rotational number J 
and the gas temperature suggested by Varanasi. (See ref. 23.) 

In this work the cross sections given in reference 24 are employed. These cross 
sections allow for the dependence on J but not on temperature . Allowing for the tem- 
perature dependence suggested in reference 23 reduces the calculated gain coefficient by 
25 percent. Thus, the discrepancy does not result from failing to use the proper gain 
coefficient or the proper cross section. 

The discrepancy can be attributed to a number of factors, the most important of 
which is that the theory assumed a rather clean experiment in which impurities such as 
H 2 O, CO 2 , and the nitrogen oxides were nonexistent. This, of course, is not a realistic 
assumption. Another reason for the discrepancy is that O 2 was treated as an inert gas 
in the theory; O 2 was added in the experiment to prevent the dissociation of CO and the 
formation of carbon. Because the calculations employed the experimentally determined 
electric field and current, the effect of O 2 on performance cannot be attributed to a change 
in the electron distribution function resulting from changes in E/N which accompany 
the addition of O 2 . (See ref. 25.) It is believed that the presence of O 2 results in the 
formation of NO and this can deactivate the excited levels of CO. 

In addition to these effects, recent analysis of gain measurements has indicated that 
measured values of gain coefficient can be substantially below the true value ; this condi- 
tion takes place in the presence of internal gain saturation and during the employment of 
large collector apertures. (See ref. 26.) The situation was further compounded by the 
fact that for low CO concentrations, the laser operation was not very stable. This insta- 
bility could be a result of an ionization instability. (See refs. 27 and 28.) 

CONCLUDING REMARKS 

The theory presented here suggests that low transitions in CO laser systems can be 
achieved at low flow rates (or pressures) for rather low CO partial pressures relative to 
those of He and N 2 . The primary role of N 2 appears to be the reduction of the electron 
temperature to a level where the transfer of electric power to the lower vibrational states 
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is optimum. The trends predicted by the theory are substantiated by experiment. How- 
ever, further analysis is needed to assess the role of impurities that inherently exist in 
any given experiment. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
March 5, 1976 
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APPENDIX A 


VIBRATIONAL RATE COEFFICIENTS 

The rate coefficients for vibration-translation (V-T) and. vibration-vibration (V-V) 
exchange collisions used in this analysis are calculated as the product of a transition 
probability (TP) and an elastic collision frequency. The expressions contain variable 
parameters which may be adjusted to obtain the best possible agreement with experi- 
mental rate measurements. 

The V-T rate expression used is based on the Schwartz -Slawsky-Herzf eld (SSH) 
theory (ref. 29, pp. 260-350) with suitable modifications. A typical V-T exchange colli- 
sion can be described by 

AB(v = n) + C 5= AB(v = n - 1) + C + AE(n) (Al) 

where A, B, and C represent different molecules and AE(n) is the energy lost in 
the collision which is converted into heat. In SSH theory, the vibrating molecules are 
treated as harmonic oscillators. Keck and Carrier (ref. 30) substituted a Morse oscil- 
lator model and also introduced an empirical fit for the adiabaticity factor (which shall be 
defined later). This latter modification allows the modeling of collision types ranging 
from impulsive to adiabatic energy exchange. Bray (ref. 31) was the first to incorporate 
these modifications of SSH theory into an analytical study of vibrational relaxation pro- 
cesses. Since that time, many theoretical CO laser studies (refs. 4 to 6 and 13) have 
made use of his work. 

A further modification to SSH theory was made necessary by a consideration of 
intermolecular interaction potentials between colliding molecules. Simple SSH theory 
assumes only an exponential repulsive potential between the molecules, but, in reality, 
an attractive force is also present. This attractive well can have a significant effect 
especially at low temperatures. (See refs. 32 and 33.) Shin (ref. 32) has shown that the 
TP for a Morse potential, which has an attractive well and also provides a fairly good 
approximation to the true intermolecular potential, can be obtained from the TP of the 
exponential repulsive potential of SSH theory by means of a multiplicative correction 
factor. 

After applying these corrections to the expressions for the V-T rates from SSH 
theory, the following results are obtained: 
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The adiabaticity factor expressed as 

F(y)=|exp(-|.)|3-exp(-|^)] 
for 0 = y = 21.622 and as 
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(A3) 


(A4) 


F(y) = exp(-3y2/3) 


(AS) 


for y > 21.622. Equation (A4) is Keck and Carrier’s (ref. 30) empirical fit and equa- 
tion (AS) comes from original SSH theory. The exponential term involving T*^ is Shin's 
(ref. 32) correction factor. The well depth €g^ can be expressed as 


2^ 


(A6) 
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The collision frequency Zg^ is calculated from the Lennard-Jones potential parameters 
of the colliding molecules (ref. 34) 


where 


^st 


- 4 " 


( 2 , 2 ) 



SffkT 

Mst 


(A7) 


‘^st 



(A8) 


The combining rates found in equation (A6) are from Good and Hope (refs, 35 and 36). 
The term n^^’^^T) is given by an empirical curve fit (ref. 37) 


n(2>2)(T) = + 0.52487 exp(-0.7732T) + 2.16178 exp(-2.43787T) 

- 6.435 X 10"'^ sin . 7.2737lj 


(A9) 


The parameters appearing in equations (A6) and (A8) are given in table V, 


TABLE V.- LENNARD-JONES POTENTIAL PARAMETERS 
[From refs. 34 and 3^ 




Is' 

cm 

K 

eV 

3.590 X 10-8 

110.0 

14.013 

3.617 

97.0 

15.576 

2.576 

10.22 

24.586 


In SSH theory, the energy -dependent TP's are obtained from a quantum mechanical 
distorted wave approximation (refs. 29 and 33). The temperature -dependent TP's are 
obtained by averaging the energy -dependent TP's over a thermal Boltzmann distribution. 
At the cryogenic temperatures at which most lasers operate, the energy -dependent TP's 
are valid. However, the thermal averaging for the temperature -dependent TP's assumes 
that the energy lost by the oscillating molecules is much less than the kinetic energy of 
the colliding molecules. At low temperatures, this assumption is not valid and the ther- 
mal averaging should be done over the Boltzmann distribution with numerical integration 
techniques. (See refs. 33 and 39.) However, this procedure is too costly to be used in 
CO laser studies. 
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Experimental measurements of the V-T rate coefficient for CO deactivation 
by He and H 2 indicate that when these coefficients are expressed in the proper vmits 
(sec“^-torr"^), they are constant below a certain temperature (ref. 40). However, if 
equation (A2) is used to plot theoretical V-T rate coefficients against temperature (after 
converting to the proper units ^), the resulting curve decreases smoothly with decreasing 
temperature until a certain critical temperature T®* is reached and then begins to 
increase with decreasing temperature. Therefore, equation (A2) is used to compute the 
rate coefficients above the critical temperature. Below that temperature, the rates are 
calculated from 


VTf (T) = 


,st 
crit 


VT 


St/rpSt \ 

crit/ 


(AlO) 


rSt 


It has been found that for a given set of constants Kyrp and ^st. the critical tempera- 
ture decreases slightly as the vibrational level quantum number increases. In this 

st 

analysis, only the T°^.|. value for v = 1 is used; this procedure has the effect of 
slightly overestimating the V-T rates for the higher levels at temperatures below T®* 

st st 

The constants Kyrp and and the values are presented in table VI. 

They were obtained by fitting equations (A2) and (AlO) to the available experimental rate 
data. (See refs. 40 to 47.) 


TABLE VI.- EMPIRICAL CONSTANTS FOR V-T RATE COEFFICIENTS 



i^st 

ft^vT 

^st’ 

0 

A 

rpSt 

^ crit’ 
K 


CO-CO 

4.4516 X 10-5 

0.203341 

135 


CO-He 

7. 1544 

.29728 

152 


N 2 -N 2 

2.8557 

.20328 

134 


N2-He 

15.476 

.3165 

145 



Figures 10 to 13 show the theoretical rates for CO deactivated by CO and He and 
for N 2 deactivated by N 2 and He together with the data points used in the curve -fitting 


^Equation (A2) has units of cm^/sec and can be converted to sec ^-torr ^ by 

iSt / 


^ 1 nocc r 


1.0355 X 10' 


28 


APPENDIX A 


process. The solid circle points in figures 10 to 12 do not represent actual data points 
but were generated from the following empirical curve fits: 


prco.co<»‘">-=»-=) = 1 o(89'3T-1''3-9.64) 
for 1372 K g T S 2915 K (ref. 41), 


(All) 


pT^Q jjg(atm-sec) = exp(87T"^^^ - 19. l) 


for 580 K s T S 1500 K (ref. 42), and 

= lo(^®^'^ ^/^-11.24) 


(A12) 


(A13) 


for 1900 K 2 T S 5600 K (ref. 43). These rate expressions are in the form of a relaxa- 
tion time pTg^ which is the standard form for giving experimentally measured results. 
This relaxation time can be expressed in terms of the deactivation rate of the first level 
of a diatomic molecule (in cm^/ sec) as 


VT 


St 

1 


1.36236 X 10"22 t 


P^st 


exp 



(A14) 


The curve fits given by equations (All), (A12), and (A13) are given in references 41 to 43 
instead of tabulated experimental results. 


No experimental measurements of the rates of CO deactivated by N 2 and N 2 deacti- 
vated by CO are available. Since CO and N 2 have the same molecular weight, SSH theory 
(ref. 29, pp. 260-350) predicts that the parametric constants appearing in the rate expres- 
sions for these reactions will be the same as those for CO deactivated by CO and N 2 deac- 
tivated by N 2 . Thus, it is assumed that 


CO-N 9 

VT^ 


^ y^CO-CO 


(A15) 


and 


N 9 -CO 

VT ^ = VT 

V V 


N2-N2 


(A16) 


All backward V-T rate coefficients were obtained from the principle of detailed balancing 
which gives 
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B vrf = vxf exp {- (A17) 

The rate coefficients for vibration -vibration exchange collisions are based on the 
work of Jeffers and Kelley (ref. 48). They proposed that V-V TP's could be expressed as 
the sum of two contributions, one a long-range multipole interaction and the other a short- 
range repulsive force. 

Single quantum V-V exchange collisions can be described by 

AB(v = n) + CD(v = m) st AB(v = n - 1) + CD(v = m + 1) + AE(n,m) (A18) 

Energy from the vibrational mode of one diatomic molecule is transferred to the vibra- 
tional mode of another diatomic molecule. If the energy spacings of the vibrational tran- 
sitions of the two molecules are equal, the energy transfer process is exactly resonant 
and AE(n,m) is zero, that is, no energy is lost as heat. The quantity AE(n,m) is 
known as the resonance defect. It is the energy transferred to the translational mode of 
the gas when the energy spacings of the transitions in the two molecules are not equal. 

The V-V TP for exactly resonant collisions is high but decreases as the resonance defect 
increases. V-V transitions with small resonance defects are dominated by the long-range 
multipole interactions of the colliding molecules (dipole -dipole, dipole -quadrupole, etc.). 
(See refs. 48 and 49.) The short-range repulsive portion of the intermolecular potential 
dominates transitions far from resonance. Thus, Jeffers and Kelley (ref. 48) suggested 
that the V-V TP's covdd be expressed as the sum of these two contributions, an idea that 
was independently developed by Caledonia and Center (ref. 50). They utilized the Sharma- 
Brau theory (ref. 49) to calculate the long-range dipole-dipole contribution to the CO-CO 
V-V TP's and determined the short-range repulsive force contribution by numerically 
averaging the energy -dependent Rapp-Englander-Giolden transition probabilities (ref. 51) 
over a thermal Boltzmann distribution. When this theory was applied to CO-CO V-V 
transitions, good agreement with experimental data (ref. 52) was obtained for tempera- 
tures around 300 K. The validity of the basic premise of this model has been questioned 
(refs. 53 and 54), but a simplified version of it is widely used because it agrees very well 
with available experimental data. 

The V-V rate model used in this analysis is a simplified version of that developed 
by Jeffers and Kelley. It was first used by Rockwood et al. (ref. 55) who showed that the 
extremely complicated rate expressions of Sharma and Brau can be approximated very 
well by a simple Gaussian function of the resonance defect. This function is character- 
ized by two parameters - its magnitude at exact resonance and its half -width. Rockwood 
et al. (ref. 55) originally assumed constant values for these two parameters, and hence 
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their expression was valid only at a single temperature. However, the proper tempera- 
ture dependence (ref. 49) of both the magnitude and half -width of the Gaussian function 
(ref. 56) is easily introduced. Rockwood et al. (ref. 55) also calculated the short-range 
repulsive force contribution using a simple expression from modified SSH theory for 
V-V transitions (refs. 30 and 31). 

The V-V rate coefficient expressions, as used in this analysis, are 
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(A22) 


and bg^, 6g, 6^, T*^, and F(y) may be obtained from equations (A3) to (A5). The 

vibrational matrix element of the dominant electric multipole moment of molecule s 

is R® ^ i/l^l 0 • dipole moment is dominant for CO and the appropriate dipole 

matrix elements may be found in table I. The electric quadrupole moment is dominant 
for N 2 and its appropriate matrix elements can be calculated from 
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Ng 
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(A23) 


where 77 = This expression is for g. linear quadrupole moment with a Morse oscil- 

lator (ref. 57) and it agrees very well with the resxilts of Truhlar (ref. 58). 

st 

The variable parameters A^^, and were determined by curve -fitting 

equation (A19) to experimentally measured rate data for CO-CO and CO-N 2 V-V transi- 
tions. (See refs. 46, 52, and 59 to 66 .) The values for these constants and for ^st are 
given in table VII. 


TABLE VII.- EMPIRICAL CONSTANTS FOR V-V RATE COEFFICIENTS 



■^st 

Wst 

K-st 

K-VV 

0 

A 

CO-CO 

CO-N 2 

N 2 -N 2 

1.16681 

1.42188 

64.155 

47.6964 

3.06156 X 10-'^ 
6.954 X 10-8 
3.06156 X lO-"^ 

0.203341 

.20331 

.20328 




No experimental data were available for N 2 -N 2 V-V transitions, but since CO 
and N 2 have the same molecular weight and very nearly equal vibrational energy level 
spacings, SSH theory (ref. 29) predicts that 


t^CO-CO 

Kw 


« K. 


N2-N2 


W 


(A24) 


However, when this assumption was used, false inversions in the number densities for N 2 
were obtained on the upper vibrational levels even when relatively large numbers of vibra- 
tional levels were carried in the calculations. Thus, this constant was decreased by a 
factor of 10. The values of fcO-CO ^N 2 -N 2 taken from the V-T rate expres- 

sions and 


^C0-N2 = \/^C0-C0%2-N2 

This expression is based on the combining rule for Lennard-Jones potentials given in 
equation (A 8 ) and simple SSH theory (ref. 29). The long-range quadrupole -quadrupole 
contribution to the N 2 -N 2 V-V TP's was found to be negligible and was not used. 
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The V-V rate coefficients given by equation (A19) automatically satisfy detailed 
balancing, that is, 




ss 


.ss 


and 


VV®^ - - = W®^ 

''''v*+l,v-l **v,v' 


exp 



(A26) 


(A27) 


Figures 14 to 18 show the theoretical W rates for CO-CO, CO-N 2 , and N 2 -N 2 tran- 
sitions together with the experimental data points used in obtaining the cvirve fits. In 
figure 14 the data of Wittig and Smith (ref. 59) was not used in the curve fitting because 
their experimental techniques have been questioned (ref. 60). In figure 17 the data from 
Sato et al. (ref. 61) and von Rosenberg et al. (ref. 62) were not used because of the large 
scatter in their results. 
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ELECTRON DISTRIBUTION FUNCTION AND EXCITATION RATES 

To obtain the electron impact vibrational excitation and ionization rates which 
appear in the governing equations, an electron distribution function is obtained from a 
solution of the electron Boltzmann equation (refs. 28 and 67) and these rates are then 
calculated in terms of that fxmction. 

The electron Boltzmann equation used in this analysis is an adaptation of the one 
for electrons in a pure gas (ref, 68) and is derived under the assumption that the density 
of the gas is high enough to insure that the electron density gradients may be neglected 
(ref. 68). It can be written as 



This is a difference -differential equation which requires a numerical solution tech- 
nique. To facilitate the solution, Phelps and his associates (refs. 68 and 69) have devel- 
oped approximations for the large number of inelastic rotational terms in this equation 
which are valid when the average electron energy is much greater than the gas tempera- 
ture. This simplification is extremely important since, in order to obtain an accurate 
solution of equation (Bl), the energy difference between mesh points must be taken 
smaller than any of the Ugj. For CO, the extremely large number of mesh points that 
would be required would make a solution of the equation impractical. When rotational 
excitation is dominated by the electric dipole moment as in CO, the rotational inelastic 
terms can be approximated by (ref. 68) 

Ns^sD «“)] (B2) 
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where 

_ 8ir „2„2 
<^sD - "3" ^O^sD 


(B3) 


For N 2 , where these processes are dominated by the electric quadrupole moment, the 
rotational inelastic terms are approximated by 


NsKsQ !;[» «“)] 


where 


^sQ = 4Bgo(TsQ 

_ 8tt „ „2 
*^sQ “ TS ^O^sQ 


(B4) 


(B5) 


In addition, it is assumed that the populations of excited states are much smaller than 
that of the ground state. Thus, the last two terms of equation (Bl) can be neglected. 
With all these approximations, equation (Bl) becomes 


A|^A(u)|t + B(u) fj - C(u) f + D(u) = 0 

where 

B(u) = XQ 2 ^ ^ 

> 

C(u) = u^^s^QsJ^’^) 

s J 

^ ^s I + "s j) + «sj) QsjI'^ + ^sj) 

s J J 


(B6) 


(B7) 
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and 


N 


Ql " ^ ^s^sm 
s 

M = ^ 
s 

Q2 = ^ /3sQsm 
s 

Kd = 2 ^S^SD 
s 

Kq = 2 /Ss^sQ 


M 

nic 


X = 


s 

2n^ 

M 


Uip = 


e 


(B8) 


All the momentum, vibrational, and electronic excitation cross sections and the ionization 
cross sections which appear in equation (B6), except those for the electronic excitation 
of CO, can be found in tabulated form in reference 69. The CO electronic cross sections 
are calculated from the semiempirical formulas of Sawada et al. (ref. 70). 

To facilitate the integration of equation (B6), it can be written as two first-order 


differential equations 


df _ g(u) - B(u) f(u) 
du A(u) 

(B9) 

^=C(u) f(u) -D(u) 
du 

(BIO) 

where 


g(u) = A(u) ^ + B(u) f(u) 

(Bll) 
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The initial conditions on f and g are 
^(%iax) “ ® 


(B12) 


where Iq is a very small constant (-10“^®). The value of is chosen large 

enough to insure that f(u) can be assumed to be essentially equal to zero for u > Umax* 
Equations (B9) and (BIO) are integrated backward from u = Uj^g^ to u = 0 by using a 
constant-step-size, fourth-order Adams-Moulton predictor -corrector method. The step- 
size value is dependent on the E/N value. The solution for f(u) is normalized so that 



1 


(B13) 


As a check on the accuracy of the solution for f(u), the electron energy equation is 

/ Op \l/2 

calculated as well. This equation is obtained by multiplying equation (B6) by -^1 u du 

\me/ 

and integrating it over all energies (ref. 68) and it is a statement of the fact that the power 
input to the electrons from the electric field must be exactly balanced by the power losses 
of the electrons from elastic and inelastic collisions. It can be written 




u Qsj(u) f(u) du 


(B14) 


where 


w = (iuN)^ 


(B15) 


and 


^N=-ij2e r - fdu 


(B16) 
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If the left and right sides of equation (B14) differ by more than 1 percent, the integration 
for f(u) is repeated with a smaller step size. The electron temperature is given by 

(B17) 

Thus, for a given E/N and a given composition, an electron temperature can be deter- 
mined from the distribution function solution. By obtaining solutions for several E/N 
values for a given composition, a table of electron temperatvires for various E/N values 
can be constructed. Interpolation in this table gives the proper electron temperature 
value for the conditions in the laser cavity. 

The current density can be related to the electron drift velocity by 


j = eNgW 

(B18) 

or 


j = eNgMN| 

(B19) 

Equation (B19) can be used to calculate the electron number density for the given set of 
conditions under consideration. 

The rate for electron impact excitation from the vibrational ground state can be 
written as 


(B20) 


The rate for electron impact deexcitation from vibrational state v to the ground state 
is 


‘’vO = lo” ““ 

where the principle of microscopic reversibility is used to obtain Qg^ from 
'iQsv^") = (" + "sv)Qsv(" + "sv) 


(B21) 


(B22) 


Results of distribution function calculations for two gas mixtures are presented in 
figures 19 to 21. Figure 19 shows the variation of Pq^ (for v = 1, 4, and 8) with E/N 
for the case 
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— 1,96 X 10 g/s6C 
rhgg = 3.91 X 10"'^ g/sec 
and figure 20 shows the same results for the case 
^QQ = 1*92 X 10"^ g/sec 
rhjg ^2 = 3.81 X 10“® g/sec 
rhjjg = 3.83 X 10 g/sec 


Figure 21 shows the variation of electron temperature with E/N for the 
The constants which appear in equations (B2) to (B5) are 


®CO,0 = 2.4x10-4 ev 
0Co,D = 4.6 X 10-2 
®N2,0 = 2.5 X 10 “4 eV 

%2,Q= 1-04 ea2(67) 
aQ = 5.29 X 10”^ cm 


= 13.6 eV 


(B23) 


(B24) 


same two cases. 


(B25) 


39 


APPENDIX C 


GAIN COEFFICIENT 


The expression used in this analysis for the gain coefficient for a given P(J) branch 
transition where the rotational quantum number changes from J - 1 to J as the vibra- 
tional quantum number changes from v to v - 1 is due to Lacina (ref. 7). It is given 

by 


r. 


V,J 


_i JS 

STTi^fkT 


R. 


CO 

V,v-1 


R 


CO 

1,0 


^CO,v 


exp 


_ b£- B J( J - 1) 
kT V ^ ^ 


- Nco,v-l®v-l 




(Cl) 


where 




(C2) 


The expression for G^q(v) is given by equations (16) and the values of A^^^ and 


pCO /pCO 
\v-l/^l,0 


are given in table I. 


The derivation of equation (Cl) assumes that the rotational levels are in equilibrium 
at the given gas temperature. The effects of both Doppler and pressure broadening (ref. 7) 
are included in the expression for the line shape factor 


where 


S = — — exp(x^) erfc(x) 



A.V 


c 


b. 

2tt 


and (ref. 68) 


f 


c 



CO,s 
s^opt 


N„a 


/ 2kT \^/2 
\^CO,sj 


(C3) 

(C4a) 

(C4b) 

(C4c) 

(C4d) 
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The optical cross sections are taken from Williams et al. (ref. 24) 


CO-CO 

%pt 

CO-No 

°"opt 


= 8.0 X 10“^® cm^ 
= 7.35 X 10“^^ cm^ 






CO-He 

°^opt 


= 3.0 X 


10“^® cm^ 


j 


(C5) 


These cross sections decrease with increasing rotational quantum number J. Equa- 
tions (C5) give values corresponding to J = 11. 

The frequency for a given P(J) transition is given by 

■-v.J = + '’vW - 1) - fv-lWi] (C«) 

where 

Fy(J) = ByJ(J + 1) - Dyj2(j + 1)2 (C7) 

and 

(C8) 
(C9) 

The units of F^, B^, and are cm"^. The constants in equations (C8) and (C9) are 

taken from Patel. (See ref. 15.) 


By = 1.93141 - 0.01752(v 


+ 2.96 X 10' 


1\2 




Dy = 6.18 X 10"® - 1.76 X 10"^ ('"+!) 
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Case 7: = 26.78 x 10"^ g/sec 
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Case 4: % = 3.81 x 10"^ g/sec 
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Figure 2.- Influence of N 2 fraction on maximum small 
signal gain coefficient. 
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Figure 3.- Influence of N 2 fraction on relative lasing intensities. 
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Upper vibrational level 


Figure 5. - Influence of CO fraction on relative lasing intensities. 
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Upper vibrational level 

Figure 6.- Effect of pressure on maximum small 
signal gain coefficient. 
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Figure 7. - Effect of pressure on relative lasing intensities. 
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Upper vibrational level 

Figure 8.- Comparison of calculated and measured relative intensities. 
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Figure 20. - CO electron vibrational excitation rate coefficients as 
function of E/N for CO-N 2 -He mixture of case 4. 
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Figure 21.- Electron temperature as function of E/N for 
CO-He mixture of case 1 and CO-N 2 -He mixture of 
case 4. 
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NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 

Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. Also includes conference 
proceedings with either limited or unlimited 
distribution. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace, 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 
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